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We apply the transfer matrix DMRG to frustrated quan- 
tum spin chains, going down to T — 0.025 while being in the 
thermodynamic limit. The incommensurability problem of 
exact diagonalization and the negative sign problem of quan- 
tum Monte Carlo vanish completely. To illustrate the method, 
we give results for chains with next-nearest neighbor frustra- 
tion and the delta chain, which has been a testbed for many 
thermodynamic methods. By comparison, the DMRG proves 
to be an extremely powerful method for the old problem of 
the thermodynamics of frustrated systems. 
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In the recent past, there has been enormous interest in 
one-dimensional quantum systems, concentrating on the 
Hubbard and Heisenberg models. As analytically exact 
methods such as the Bethe or the matrix product ansatz 
can capture only some of the relevant models, particular 
attention has been paid to the development of numerical 
methods. 

Since one-dimensional strongly correlated systems 
show strong quantum and thermal fluctuations, the most 
interesting physics is to be expected close to or at T = 0, 
where only quantum fluctuations survive. To investigate 
such systems in the low temperature regime, several nu- 
merical techniques are used, mainly exact diagonaliza- 
tions, quantum Monte Carlo, and the DMRG algorithm 

0- 

In many cases, these methods yield completely sat- 
isfying, mutually complementary or supporting results. 
Some systems, however, have evaded thorough analysis 
so far. One of the strongest numerical challenges comes 
from frustrated quantum spin models 0] . These models 
are both of theoretical and experimental interest, which 
at the moment concentrates on the famous CuGe03 ||. 




FIG. 1. Hamiltonian studied by the transfer matrix 
DMRG. The dashed line indicates spins grouped into a new 
object. 

In this letter we show how the transfer matrix DMRG 
can be applied to frustrated quantum spin systems, such 



that thermodynamic quantities can be very reliably ob- 
tained down to T sa 0.01 for infinite system sizes, i.e. in 
the thermodynamic limit. There is no principal restric- 
tion of the systems that can be considered. To discuss 
how this overcomes the difficulties of other methods, we 
present numerical results for several physically relevant 
frustrated quantum spin systems. 

Let us first recall the other standard techniques. Exact 
diagonalization techniques either calculate the full spec- 
trum of a very small system or extract information from 
the low-lying states of somewhat longer systems (J] . In 
both cases, finite size effects are present, aggravated in 
frustrated systems by incommensurability effects. 

Both the quantum transfer matrix technique J5| and 
quantum Monte Carlo start from a Trotter decomposi- 
tion which maps the one-dimensional quantum problem 
to a two-dimensional classical problem. After a chequer- 
board decomposition H = H 1 + H 2 ||, the partition 
function is given by 

Z m = Tr [ e -PBi/™ e -l*B a /myn } (!) 

which becomes exact for the Trotter number m — > oo. In 
the quantum transfer matrix method, a transfer matrix 
relating states on sites i and i + 1 is written down ex- 
actly and thermodynamic quantities are calculated from 
its maximum eigenvalue. Matrix size grows exponentially 
with to, restricting to a rather coarse Trotter decompo- 
sition, making access to low temperatures difficult; how- 
ever, results have no finite size effects. Quantum Monte 
Carlo can access much larger Trotter numbers, while 
working on finite system sizes. In frustrated spin sys- 
tems, the essential limiting factor comes from the well- 
known negative-sign problem: some of the matrix ele- 
ments of the infinitesimal evolution operator acquire neg- 
ative signs, but must be interpreted as positive probabil- 
ities. Only partial elimination of this problem has been 
achieved 

Recently, the application of the DMRG to both classi- 
cal [[nj and quantum transfer matrices has been pro- 
posed. For the Heisenberg chain, thermodynamic quan- 
tities can be calculated very well in fact, the re- 
sults are also superior to those obtained from related at- 
tempts to apply DMRG techniques to thermodynamics 
p5[ . However, for the systems studied, other methods 
had already yielded very good results ]T^ |. 

Here, the transfer matrix DMRG is applied to the more 
challenging frustrated spin-^ Hamiltonian 
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•h^2i ' &2i+2 + ^Sai-l ' S2i+1 + iJS'f, (2) 

as shown in Figure |l|. It comprises various frustrated spin 
models of interest. 

The modifications to the original transfer matrix 
DMRG for an unfrustrated antiferromagnetic chain 
|l3| , fl4f are very simple in principle. First, we rewrite 
the frustrated spin chain as a spin ladder problem. If 
we group two neighboring spins- ^ into one object with 
4 states, every second nearest neighbor interaction is 
absorbed; the others and the next-nearest neighbor in- 
teractions become effective nearest-neighbor interactions. 
Similar in spirit to the restructuring method of Munchisa 
and Munehisa flCf| and Nakamura's procedure for allevi- 
ating the negative sign problem in quantum Monte Carlo 
, we are now free to choose any basis of the two-spin 
object which consists of eigenstates of S z , to maintain the 
good quantum number (l7) of transfer matrix DMRG, 

£(-l)W, (3) 
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where the sum runs along the Trotter direction. 

For a S = 1/2 problem, these are essentially the Ising 
basis | |T)i | TJ-)j I IT)) I 14) an d the singlet-triplet basis 
I TT>, (I U) + UT))/V2, (| U)-UT»/\/2, | ||). Numer- 
ically, we find that both yield the same results, though 
the latter has no update of its internal bond in a transfer 
matrix step and might therefore be expected to be more 
precise. All results given below were obtained with the 
singlet-triplet basis. 

The transfer matrix DMRG is governed by the param- 
eters M, the number of states kept, and (3q, the initial 
temperature and at the same time the stepwidth of the 
values of (3 considered during the DMRG iterations. As 
the two-spin object has 4 instead of 2 states, for given 
computational resources, M has to be halved. To esti- 
mate the effect on DMRG precision (and also check the 
correctness of the program), we applied both the origi- 
nal and our transfer matrix DMRG to the unfrustrated 
Heisenberg chain. We find that the new version shows 
much better results for M/2 than the original one for 
M states kept, as every second bond is taken into ac- 
count exactly. This greatly enhances the applicability of 
this method also to frustrated chains. To give one result, 
the internal energy per site of the isotropic Heisenberg 
model should go to -0.4431 for T — > 0. For M = 32, the 
original transfer matrix DMRG extrapolates —0.448 for 
(3q = 0.2, ours gives -0.445 for M = 16 only. 

There are several advantages of the transfer matrix 
DMRG for frustrated spin systems. As opposed to the 
exact diagonalisation of small systems, we have no in- 
commensurability problem, as we work in the thermody- 
namic limit right away. Unlike quantum Monte Carlo, 



there is no such thing as a negative sign problem; also 
quantum MC results are obtained in finite systems, al- 
beit often rather large ones. Compared to the recent 
progress made by Nakamura ||, DMRG works for all 
frustrated spin systems; furthermore, there is no problem 
with the susceptibility. Quantum transfer matrix meth- 
ods are closest in spirit to the transfer matrix DMRG. 
They relate to each other exactly as T = DMRG re- 
lates to exact diagonalisation. Like T — DMRG that 
can treat much longer systems than exact diagonaliza- 
tion, here much higher Trotter numbers can be obtained: 
Typically, in quantum transfer matrix problems one is re- 
stricted to less than 10 Trotter slices when dealing with 
frustrated problems. The transfer matrix DMRG has no 
fixed number of Trotter slices; it increases with (3. Con- 
sider T = 0.025, or (3 = 40 as a typical low temperature 
one might wish to reach. For a realistic starting temper- 
ature 0o — 0.05, this implies 800 Trotter slices at that 
temperature, or 800 DMRG steps. While maintaining the 
advantages of the quantum transfer matrix, the problem 
of convergence of results in the Trotter number can be 
virtually eliminated, as the results presented here show. 
However, the problem of convergence in the number M 
of states kept has also to be considered. Working with 
M < 32, which is still possible on a workstation, we find 
that down to T m 0.05 results are very well converged; 
for lower T the trend is very clear (see as an example the 
inset of Figure ||). 

In order to test the transfer matrix DMRG, we have 
first applied it to the so-called delta chain, because it 
exhibits a rather special low-temperature behavior [ p0[ , 
of which some approximative analytic results are known 
pjj . Numerical methods |?],|],^0| have been applied to 
this chain with varying success, which makes this a good 
starting point for a comparison of results. 

The delta chain consists of spin triangles, which are 
linked through the triangle basis spins (Ji = Ji = J3 = 1, 
J4 = 0) . It can be understood as a frustrated chain with 
next-nearest neighbor interactions, of which every second 
has been taken out. It has the same doubly degenerate 
ground state as the Majumdar- Ghosh model, with spin 
singlets sitting cither on all left or on all right sides of 
the triangles [Q . Numerical studies reveal that the delta 
chain gap A w 0.22 is much smaller than the Majumdar- 
Ghosh model p9| gap, which suggests a rather delocal- 
ized excitation |20|| . On the other hand, exact diagonali- 
sation shows an almost flat dispersion relation, indicating 
a localized excitation; these results are however strongly 
size-dependent [^0| . It has been suggested that there are 
two types of excitations in the chain, kinks and antikinks 
pl|,p2[ . Elementary excitations can be imagined as do- 
main walls between ground state segments with a singlet 
on the left to singlet on the right or right to left flip. In 
the former case, it is energetically favorable to insert one 
triangle with a singlet on the basis. This is low in energy, 
being an eigenstate of the local triangle Hamiltonian, but 
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localized. In the latter case, one triangle has no singlet 
at all, which is high in energy. However, this triangle 
can be delocalized. Obviously, kink and antikink excita- 
tions must alternate. Thus the spectrum finds a natural 
explanation. 

For the specific heat, the presence of two different ex- 
citation types with different energies implies a two-peak 
structure, with a conventional high-temperature peak 
and an additional low-temperature peak. 

We have calculated the specific heat C V (T) by numer- 
ical differentiation of the internal energy which was di- 
rectly evaluated by the DMRG, starting from tempera- 
tures (3o = 0.5, 0.4, 0.2, 0.1 and 0.05, while retaining 
M = 16, 24 and 32 states, down to temperature /3 = 40. 
This implies maximum Trotter numbers m = 80, 100, 
200, 400 and 800. 
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FIG. 2. Specific heat of the delta chain for 32 states kept 
and various /3o- Note the different convergence behavior of 
the high- and the low-T peak. 




FIG. 3. Low temperature peak of the specific heat of the 
delta chain (M = 32). Curves for /3q — 0.1 and 0.05 are 
almost indistinguishable. 



We find that the low-temperature peak in the specific 
heat is converged both in /3q and M; we extract ther- 
modynamic limit values of T pea fc = 0.082(2) for the peak 
temperature and a value C™ ax — 0.184(1). Calculating 
the specific heat from a variational ansatz, Nakamura 
and Kubo |H] predict this second peak at T s» 0.05, 
with a height of C v ss 0.24. Numerically, Kubo finds 
(jmax gjg at x peak w 0.12 from a quantum transfer 
method pOfl ; these are unextrapolated results from his 
maximum Trotter number. Otsuka finds C™ ax « 0.20 at 
Tpeak ~ 0.12 from a recursion method |23|. Monte Carlo 



studies have reached the high-temperature slope of the 
low-T peak of C v , but could not localize and estimate 
it reliably ||. The DMRG therefore clearly gives the 
quantitatively most accurate results; Nakamura's modi- 
fied Monte Carlo || should in principle be capable of a 
similar precision for the problem under study; it must be 
stressed that this is because the delta chain fits his neg- 
ative sign problem elimination formula which is however 
not universal. The DMRG exploits no special feature of 
the studied frustrated system and is thus more versatile. 

Let us comment on the convergence of the specific heat 
with 0o — > 0. The underlying internal energy is uniformly 
underestimated for finite /?o, converging as /?q. However, 
as can be seen in Figure ||, for coarse Trotter slicing, 
the high-temperature peak is overestimated, while the 
low-temperature peak is underestimated (a similar phe- 
nomenon is also visible in Ref. [p0[). This can be under- 
stood in terms of a "sum rule" . Obviously, 



s(T = oo) - s(T = 0) = / 
Jo 



C V (T) 
T 



dT. 



(4) 



Now all the classical models resulting from the Trotter 
decomposition have vanishing zero-temperature entropy, 
while they all map free spins for T — > oo, i.e. s(T — oo) = 
In 2, yielding a constant value of the integral for all (3q. 
Thus the errors due to the Trotter decomposition have to 
compensate each other; a coarse decomposition will miss 
low-T features and overemphasise high-T features. The 
shoulders for (3 = 0.5 and 0.4 in Fig. || we interpret as a 
low-T cutoff of C v due to large (3$. 

For a comparison, the specific heat (Fig. ||) of a typical 
next-nearest neighbor frustrated chain (the spin-^ chain 
at its transition point J\ — J2 = 1, J3 = J4 = 0.2411) 
has a one-peak structure, which is converged down to 
T w 0.04 and can be linearly extrapolated to C v = 
at T = 0. For very low temperatures, the results devi- 
ate towards higher values. We find that both a too coarse 
Trotter decomposition and a too small M lead to an over- 
estimation of C v ; comparing the curves for /3q — 0.1 and 
(3q = 0.05 indicates that the remaining error is almost 
entirely due to the finite number M of DMRG states. As 
the inset in Fig. ^ shows, even for very low T the results 
converge fast with M. 

Let us finish by showing the susceptibility of the delta 
chain in comparison to the exactly known susceptibility 
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of the isotropic S = \ Heisenberg chain (Fig. ||) , both cal- 
culated by numerical differentiation of the explicitly eval- 
uated magnetisation. Whereas for the latter one finds a 
finite % for T = 0, here the susceptibility is exponentially 
activated due to the gap in the excitation spectrum. No 
particular physics is hidden there, but from the point of 
view of the transfer matrix DMRG it can be observed 
that again results converge very fast in /3q and M. 



method for obtaining the thermodynamics of frustrated 
spin or also frustrated fermionic systems. 

All calculations were carried out on a 533 Mhz alpha 
with Linux and a 16 processor T3E. 



0.3 



0.2 



0.1 




FIG. 4. Specific heat of the next-nearest neighbor frus- 
trated spin-i chain at its phase transition (M — 32). The 
inset shows the convergence with the number M of states 
kept of the f3o = 0.05 curve for very low temperatures. 
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FIG. 5. Susceptibility of the delta chain (M = 32) com- 
pared to the susceptibility of the Heisenberg chain (M — 32). 

To conclude, we have demonstrated how the trans- 
fer matrix DMRG overcomes the problem of obtaining 
the low-temperature thermodynamics of frustrated spin 
chains, while at the same time working in the thermody- 
namic limit without serious restrictions on the Hamilto- 
nians considered. The method can easily be upgraded to 
more complicated frustrated chains, where other ad hoc 
fixtures fail. We are therefore confident that the trans- 
fer matrix DMRG will emerge as the most convenient 
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